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We calculate the critical temperature for Bose-Einstein condensation in a gas of bosonic atoms 
across a Feshbach resonance, and show how medium effects at negative scattering lengths give rise 
^>«^ , to pairs reminiscent of the ones responsible for fermionic superfluidity. We find that the formation 

■ of pairs leads to a large suppression of the critical temperature. Within the formalism developed 

' by Nozieres and Schmitt-Rink the gas appears mechanically stable throughout the entire crossover 

^S) , region, but when interactions between pairs are taken into account we show that the gas becomes 

, ^ ■ unstable close to the critical temperature. We discuss prospects of observing these effects in a 

gas of ultracold ^^^Cs atoms where recent measurements indicate that the gas may be sufBciently 
, long-lived to explore the many-body physics around a Feshbach resonance. 

l/^ ■ PACS numbers; 03.75.-b, 67.85.-d, 67.85.Jk, 67.10.Ba 
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I. INTRODUCTION. 

In a Fermi gas there exists a smooth crossover connect- 
ing two apparently distant states, namely, the Bardeen- 
Cooper-SchriefFer (BCS) state which is found for weakly 
attractive interactions and the Bose-Einstein condensed 
(BEC) state of diatomic molecules obtained for weakly 
repulsive effective interactions fl*!. This crossover be- 
tween a superfluid of loosely bound pairs of fermions and 
a condensate of tightly bound dimers takes place in the 
vicinity of a Feshbach resonance, where the effective in- 
teraction changes from attractive to repulsive through a 
divergence of the coupling constant. What lies at the 
heart of this continuous crossover is the fact that a pair 
of fermions makes up a boson which is bound either by 
two-body effects in the BEC limit or by many-body ef- 
fects in the BCS limit. It is this composite boson that 
undergoes Bose-Einstein condensation. 

While fermions can condense only in the form of pairs, 
bosons can condense solitarily. A Bose gas at low tem- 
peratures is therefore subject to a competition between 
the condensation of atoms and of pairs. Pairing in a 
bosonic gas leads either to the formation of tightly bound 
molecules which are stable even in the vacuum, or to the 
creation of loosely bound pairs stabilized by the medium. 
The latter are reminiscent of Cooper pairs arising in 
BCS superconductivity. In this article we investigate 
the properties of a Bose gas near a Feshbach resonance 
and address important issues regarding stability of such 
a strongly interacting Bose gas. 

The BEC-BCS crossover for fermions has been recently 
explored in a series of ground breaking experiments with 
ultracold alkaline gases 3 by exploiting so-called Fes- 
hbach resonances which enable the interaction strength 
to be tuned. By contrast, early experimental attempts 



to create strongly interacting bosonic gases have been 
plagued by large losses 13\ . This is primarily because in- 
elastic collisions provoking the decay of pairs into deeply 
bound states increase rapidly in the neighborhood of a 
resonance of the s-wave scattering length a. The de- 
cay imparts a kinetic energy to the products that is in 
general much larger than the confining potential, caus- 
ing atoms to escape from the trap. Moreover, attractive 
atomic interactions may induce a mechanical instability 
at sufficiently large densities resulting in the collapse of 
the gas [^]. Both problems do not arise with fermions 
due to the stabilizing effects of Pauli blocking Q . 

However, in recent experiments a gas of 6s molecules 
composed of bosonic ^^'^Cs atoms was observed to have 
relatively small inelastic losses 0I- This stability 
against decay offers the intriguing possibility of exper- 
imentally realizing a molecular BEC with bosonic atoms 
and of studying the BEC-BCS crossover in a bosonic gas. 

The paper is organized as follows. In Sec.|TT]we demon- 
strate that a crossover exists in the case of bosons that is 
analogous to the BEC-BCS crossover occurring in Fermi 
gases. We show that the compressibility obtained within 
the Nozieres and Schmitt-Rink (NSR) formalism Q is 
positive throughout the normal phase, indicating that the 
gas may be mechanically stable. However, attractive in- 
teractions beyond the NSR formalism between long-lived 
preformed pairs of bosons can mechanically destabilize 
the gas and the elucidation of this effect is the subject of 
Sec. Ill Bl In Sec. IIIII we discuss the experimental feasi- 
bility of observing such a crossover in ^^'^Cs, and present 
our conclusions in Sec. HVl 
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II. UNIVERSAL PHASE DIAGRAM 

In principle, a two-channel model is required to com- 
pletely describe a Feshbach resonance [9]. However, suf- 
ficiently close to the resonance, a single-channel descrip- 
tion which lends itself to the NSR formalism suffices. 
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The phase diagram describing the transition from the 
normal state of a bosonic gas with density n to the con- 
densed state of either atoms or pairs is then universal, 
in the sense that it depends only on a single interac- 
tion parameter, namely, the total scattering length a di- 
vided by the average interparticle spacing n"^/"^. In the 
limit ri^/^a 0~, atoms condense at the critical tem- 
perature for an ideal gas Ta = (27r?i^/mfcB)[n/C(3/2)]^/^, 
where n is the total atomic density. In the opposite limit 
•n}/^a 0"*" there exists a deep two-body bound state, 
and atoms bind into dimers which condense at the lower 
temperature Tm — Ta/2^/'^. Within a simple two-body 
approximation (dashed line in Fig. [1]), we find a crit- 
ical temperature for pair condensation Tc that quickly 
interpolates between these two limits. However, within 
the NSR many-body approach (solid line in Fig. [1]) we 
find that Tc is strongly suppressed on the attractive side 
of the resonance. As we will see in the following, this 
suppression is due to the population of a low-energy res- 
onant two-body state brought about by many-body ef- 
fects. This is in many ways reminiscent of what occurs 
in superconductors, where electrons in a medium pair up 
below a critical temperature in the presence of an attrac- 
tive interaction. While in the fermionic case the criti- 
cal temperature falls to zero exponentially, we find that 
here Tc reaches its limiting value Ta according to a power 
law [lO|. Within the NSR approach, the gas is found to 
be stable in the normal phase at all temperatures down 
to the critical temperature Tc- In Sec. Ill Bl we will see 
that inclusion of molecule-molecule interactions leads to 
a critical temperature that is higher than the NSR result, 
and in the region close to the resonance the gas acquires 
a negative compressibility at temperatures close to T^- 



A. NSR calculation 

The results from the NSR formalism are summarized 
in Fig. [T] which shows the critical temperatures Tc and 
T° for the Bose-Einstein condensation of paired and un- 
paired bosons, respectively. The critical temperature and 
the pressure of the gas was calculated within the grand 
canonical ensemble by introducing the thermodynamic 
potential of a gas of particles of mass m at temper- 
ature T. The pressure and the density of the gas are 
related to by P = -Vl/V and n = -{l/V)dn/dfi re- 
spectively, where V is the volume and fi the chemical 
potential. Following NSR, we write as ,8] 



T^^(O) 
r(k, iujn) ' 



(1) 



The first term is the contribution of an ideal gas of atoms 
with = S^k^/2m - fj, an d /3 = l/fceT. The second 
term (rederived in Sec. IIIB[) represents the contribution 
of paired atoms where LUn = 2Tm/hf3 are the even Mat- 
subara frequencies and the two-body T-matrix is given 




FIG. 1: Universal phase diagram for bosons as a function of 
the interaction parameter where n is the total den- 

sity of atoms. The solid line shows Tc calculated within the 
NSR formalism, while the dashed line is the two-body result. 
The dash-dotted line is the critical temperature for the 
condensation of unpaired bosons in units of the critical tem- 
perature Ta of an ideal atomic gas. The phase with a pair 
condensate (PC), which lies between Tc and Tc, extends all 
the way to —l/n^^'^a = +oo. Both an atomic condensate and 
a pair condensate (AC-I-PC) exist below T^. 
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T^^(z) 



1 — a\J —zmjW- 
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where z is the energy in the centre-of-mass frame. The 
instability towards pair condensation is signaled by the 
appearance of a pole in the T-matrix at k = and iujn = 
0, which is the Thouless criterion [l^. Neglecting many- 
body effects in the first instance, we take T(k, iujn) = 
T'^^{z) with z = z(k, iujn) = if^n + 2/i - h'^kyAm. This 
two-body result is plotted as a dashed line in Fig. [TJ 

To include many-body effects, we take T(k, iuJn) = 
T^^(k,iujn), where the many-body T-matrix is [l^ 



T^^(k,*c.„) 



1 



[T2B(z)]-i-S(k,*c.„)- 



(3) 



Here, the renormalized correlation function in the 
particle-particle channel S(k, iw„) is given by 

H(k.M.)^iy (4) 



V ^ ihWn 

q 



2-q 



Ck/2+q ' 



with the Bose factor N{x) = (e''^ - 1)"^ Within this 
formalism, the Thouless criterion for condensation yields 
the solid line in Fig. [TJ 

The strong suppression of T^ for negative scattering 
lengths can be understood from the spectral function for 
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FIG. 2: Spectral function of a pair evaluated at 1/r. 
2.18, T = Tc and |k| ~ 0. The two-body (dashed) and many- 
body (solid) calculations agree at high energy, but differ sub- 
stantially at low energy. Also, a delta-function is present at 
positive frequency just below the continuum in the many- 
body case (not shown). At |k| = this delta function corre- 
sponds to the Thouless pole at tj = 0. Inset: binding energy 
of a pair ai T = Tc- The dashed line is the ideal gas limit 
£m = —h'^/ma'^, while the solid line is the NSR result given 
by Em = 2fi. 



the pairs which is proportional to Iin[T(k, oj + iO)] and 
is plotted in Fig. [5] The many-body spectral function 
(solid line) acquires a narrow peak near zero energy in the 
continuum. For small momenta, a delta-function is also 
present just below the continuum that is associated with 
the presence of infinite-lifetime pairs (not shown). Both 
features are notably absent from the two-body spectral 
function (dashed line) [9|, il6|] , and enhance the density of 
states for pairs at low energy leading to a reduction of Tc 
compared to the two-body case. In addition, populating 
these many-body resonances reduces the atomic density 
thereby reducing the atomic condensation temperature 
T". The inset shows the binding energy at Tc of the two- 
body bound state in the presence of the medium (solid 
line) which approaches the two-body result (dashed line) 
as a ^ 0+. For a < there is no bound state in the 
vacuum, and the binding energy in the medium is that of 
the pair associated with the delta-function in the spectral 
function. 

The pressure of the gas in the normal phase is plot- 
ted in Fig. [3] as a function of density up to the critical 
density for pair condensation which is denoted by a dot 
terminating the lines. We show results for n^^^a ~ 0^, 
at resonance, and for n^^^a ~ 0~. As expected, in the 
low density limit the first curve reproduces the ideal gas 
law for atoms (P — nksT), while the last yields the cor- 
responding relation for molecules (P — nkBT/2). We 
see that the compressibility (= n^^dn/dP) calculated 
within the NSR approach is always positive, implying 
that the gas in its normal phase remains mechanically 




n Km 



FIG. 3: The pressure of the gas in its normal phase as a func- 
tion of density in the NSR formalism for ~l/n^^^a ~ —10'^ 
(dashed), (solid), and 10^ (dashed-dotted). Here, A(r) = 
y^2TTh^ /rnksT is the thermal de Broglie wavelength. The 
dotted lines are the ideal gas law for atoms and molecules. 



stable throughout the whole crossover. 

Below Tc, bosons initially condense in the form of 
pairs [ll| and this region is denoted by PC in the phase 
diagram. Interactions between pairs are neglected within 
the NSR formalism, and below Tc we therefore have an 
ideal Bose-Einstein condensate of pairs. The Thouless 
criterion is then always satisfied but the total number of 
atoms found from Eq. ([T]) is too small. This is because 
the number equation misses a nonzero pair-condensate 
contribution Uc which emerges as a result of a zero- 
momentum divergence connected to the Thouless pole 
of the many-body T-matrix. For T < Tf, it is 



iV(0) aT^«(O,0) /9T"«(O,0) 



V 



\ diujn 



(5) 



The pair condensate is characterized by the BCS order 
parameter Aq = (A(x, r)) — Vo(V'(x, t)'!/'(x, r)) ^ and 
('0(x, r)) — 0. Here A(x, r) is the pairing field, Vb is 
the bare atomic interaction and '>P{'x.,t) is the atomic 
field. Noting that the Fourier transform of the BCS 
Green's function Ga(x, r; x', r') = -(A(x, r) A* (x', r')) 
can be written in terms of the many-body T-matrix as 
GA(k, zwn) — 2hT^^ (k, iujn) [4|, we find an expression 
for the order parameter |Aop = — ^„ G'a(0, ia;„)/fi/3y 
in terms of the many-body T-matrix 



lAnP 



2hN{0) 
' V 



d 



1 



9icj„TMB(o,0) 



(6) 



Now, by comparing the divergences in Eqs. ([5]) and ([6]) 
we express the gap in terms of the condensate fraction as 



|Aop = 2nc 



d 



1 



dfi TMB(o,0) 



(7) 
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Furthermore, the atomic dispersion in the pair conden- 
sate becomes a gapped BogoHubov quasiparticle disper- 
sion hu!]^ = y^S,"^ — |Aop. Thus, condensation of atoms, 
signaled by (?A(x, t)) becoming nonzero, takes place when 
the quasiparticle dispersion vanishes at zero momentum. 
The onset of atomic condensation is hence found as the 
temperature T° at which |Ao| = — This is the dash- 
dotted line in Fig. [TJ 

It is only at the temperature < Tc that an atomic 
condensate (AC) of solitary bosons appears, and coexists 
with the condensate of pairs. The phase transition sepa- 
rating the AC and AC-I-PC regions becomes an Ising-like 
quantum phase transition when T° = [H, We 
find no point at which the three phases (normal, PC, 
AC-I-PC) meet, as remains strictly smaller than Tc at 
any value of the scattering length. 
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B. Beyond NSR 

Interactions between infinite-lifetime dimers corre- 
sponding to the positive frequency pole of the many- 
body T-matrix are neglected in the NSR theory presented 
above. In order to estimate the effect of these inter- 
actions in the normal phase, we add a mean- field shift 
'T-m?nim to the binding energy of these infinite-lifetime 
dimers, and the corresponding additional interaction en- 
ergy Fn^Tmm in the thermodynamic potential. Here, 
Tmm — 47rft^anini/2m is the molecular T-matrix, a^^i is 
the dimer-dimer s-wave scattering length and is the 
density of dimers that is the contribution of the pole in 
the T-matrix to the number equation. 

To calculate the dimer-dimer scattering length, we be- 
gin by performing a Hubbard-Stratonovich transforma- 
tion and subsequently integrate out the atomic fields to 
obtain an effective action in terms of the pairing field 

hi3 



— 1 /■ f 
S^^[A,A*] = —J^ dr ydx|A(x,T)| 



^ Tr [ln(-Go 1)] ^ ^T. ' [(^0^)"] , (8) 

n=l 



where the noninteracting Green's function is a 2 x 2 ma- 
trix in Nambu space given by 



^Go(x, r;x',r')" 



-hdr 



2m 



/i 







(9) 

and the self-energy is 

;.S(x,r;x',r')= (^^,(°^^) 5(x - x')<5(r - r'). 

(10) 

The first term in Eq. ([5]) corresponds to the thermody- 
namic potential of a noninteracting gas that is the first 



FIG. 4: The renormalization factor Z (solid line). Inset: 
the dimer-dimer scattering length amm (solid line). The re- 
spective BEC limits which are recovered as a ^ 0^ are shown 
as dashed lines. 



term in Eq. ([T]). The BCS Green's function defined in the 
previous section is obtained from the terms in Eq. ([8|) 
that are quadratic in the pairing fields. They may be 
written as 

^(^) [A, A1 = -nY^A^iK ^^n) ^^^J AiK -„). 

k,n 

(11) 

Neglecting all higher-order terms and performing the 
Gaussian functional integral yields the second term in 
Eq. 11]), where the factor T'^^{0) appears due to the nor- 
malization involved in the Hubbard-Stratonovich trans- 
formation. 

By introducing the rcnormalized pairing field 
(f){'k,iuJn) — A(k, iw„)/Vz with the renormalization 
factor 



-1 



d 



1 



diujn TMB(o, 



(12) 



where is the frequency of the molecular pole in the 
many-body T-matrix at zero momentum, this quadratic 
action can be recast into the action for an ideal Bose gas 

5(2) [cj), d)*] = 0*(k, tLOn){-ih^n + ^0</'(k, ZC^„). 

k,n 

(13) 

Lowest-order corrections to NSR due to interactions be- 
tween pairs are then obtained from the fourth-order 
term of the effective action. Writing it as S''-^-' [0, 0*] = 
(n/3I^/2)|(^o|^Tmm, where (/-o = ^o/VZ, we thereby find 
the molecule-molecule scattering length 



5 



0.8 



0.6 



0.4 

T /T 

















b 



-50 -40 -30 -20 -10 10 20 30 40 50 
-1/n a 



FIG. 5: Universal phase diagram beyond NSR for the same 
parameters as in Fig. [T] The solid line is Tc computed includ- 
ing molecular interactions. The dashed-dotted and dashed 
lines are same many-body and two-body Tc results respec- 
tively shown in Fig. [1] The shading denotes the mechanical 
instability region. 




The integrand in Eq. (fH|) is sharply peaked about a; = 
thus in the BEC limit, where Z ~ IGnh'^/m^a and 
/Lt ~ —lT?/2ma^ is large and negative, the integral can 
be neglected and the molecule-molecule scattering length 
reduces to Qm m = —4a. Note that this is the same result 
found in Ref. [l3j up to a minus sign which was erro- 
neously omitted there. Plots of a„im and Z are shown in 
Fig.Sl 

The critical temperature obtained by including mean- 
field effects due to the stable dimers is plotted as the 
solid line in Fig. [SJ and is always higher than the crit- 
ical temperature predicted within the NSR formalism 
(dashed-dotted line) owing to the presence of an attrac- 
tive interaction between the stable dimers throughout the 
crossover. The pressure of the gas for densities below the 
critical density is shown in Fig. [6] at various atomic scat- 
tering lengths (solid lines). We see that the inclusion 
of molecule-molecule interactions gives rise to a region 
of negative compressibility, indicating a mechanical in- 
stability at temperatures above the critical temperature. 
The shading in Fig. O denotes the resulting unstable re- 
gion with an upper boundary given by the temperature 
for which dn/dP = 0. A striking feature of this phase 
diagram is the presence of a mechanically stable region at 
negative scattering lengths, where the gas remains in the 
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FIG. 6: The pressure of the gas in the NSR for- 
malism (dashed lines) and with molecular interactions be- 
yond NSR included (solid lines). The curves beginning at 



PA{Ty/kBT = 0.2 and 0.4 correspond to -1/; 
and 10'^, respectively. 
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normal phase far below the ideal gas critical temperature 
due to the many-body effects discussed above. 

The effects of screening of the interatomic interaction 
may be estimated by dressing the atomic interaction in 
the many-body T-matrix with an RPA bubble sum, 



1 



[r2B(z)]-i -n(O,0)-S(k,za;„)^ 



where 



n(o,o) = --^^G(k,z^„)^ 



(15) 



(16) 



k,n 



and G(\i,iujn) = h/{ihu;n — ^k)- We have estimated 
that this correction raises the critical temperature by less 
than 15% on the BCS side. The phase diagram obtained 
within a more refined calculation which includes screen- 
ing fully self-consistently is beyond the scope of this pa- 
per but should therefore not deviate substantially from 
the one presented here. 



III. CESIUM PHASE DIAGRAM 

We now turn our attention to the specific case of 
a gas of ^■^^Cs atoms near the narrow s-wave Fesh- 
bach resonance due to a d-wave molecule located at 
Bq = 47.78 G [13]. The shallow 6s bound state where 
small inelastic losses were experimentally observed is 
present just below the atomic continuum and should 
be taken into account since it crosses the 4d state that 
causes the Feshbach resonance. Thus, the two-body T- 
matrix for the -'^^'^Cs s-wave Feshbach resonance consid- 
ered here is obtained from Eq. ^ by replacing a with 
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FIG. 7: Binding energy of the 6s and 4d bare states (dashed 
lines) giving rise to the avoided crossing between the two 
dressed states (solid lines). In the close-up we show the 
Wigner threshold behavior of the binding energy near the 
continuum. 



abg(B)[I + Afj.AB/{z - S)] H, [H, [H. Here, the differ- 
ence of magnetic moments between tfie Ad bound state 
and tlie atomic continuum is A/i = 1.15/iB [l7| . and 
o-hg{B) is the background scattering length [2^. The en- 
ergy detuning from the Feshbach resonance is defined as 
6 = Aii{B — Bq) and the width of the resonance is found 
to be AB = 0.16 G [ij- 

The energy of the poles of the cesium two-body T- 
matrix are plotted as solid lines in Fig. [7] which shows 
the avoided crossing between the Ad and 6s states. The 
binding energies of the bare states are shown as dashed 
lines. Also shown is a close-up near the Feshbach res- 
onance showing the Wigner threshold behavior of the 
binding energy that is characteristic of an s-wave res- 
onance. This mirrors the behavior of the binding energy 
in the two-body limit which is plotted as the dashed line 
in the inset of Fig. [2l The many-body T-matrix is found 
by inserting the two-body T-matrix for cesium obtained 
above into Eq. Q. 

The experimentally relevant crossover is between a 
BCS-like state in the atomic continuum and a gas of con- 
densed dimers occupying the upper branch of the avoided 
crossing. To describe theoretically this metastable con- 
figuration, we neglect the contribution to the thermody- 
namic potential that arises from the deeper-lying pole 
of the cesium two-body T-matrix corresponding to the 
lower branch of the avoided crossing. The resulting criti- 
cal temperature for two densities is plotted in Fig. [51 As 
the magnetic field is increased across the Feshbach reso- 
nance, the character of the pairs changes gradually from 
stable and deeply-bound dimers to infinite-lifetime pairs 
weakly bound by many-body effects. The density of the 
pairs decreases gradually as B is increased, reaching zero 
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FIG. 8: Critical temperature for pair condensation Tc in a 
gas of ^''■^Cs for particle densities n — lO'^^cm"^ (dash-dotted) 
and lO^^cm^'^ (solid), and the corresponding instability re- 
gions (dark and light grey, respectively). For n — 10^*cm~^ 
we have Ta =79 nK and Tm =25 nK. Note that the curve has 
no kink at B — _Bo = AS = 0.16G although this is not visi- 
ble on this scale. The dashed line is the two-body result for 
n = 10"cm-^ 



as B — Bq reaches AB, where the scattering length van- 
ishes. For B ^ Bq + AB there are no longer any pairs 
and both Tc and T° reach their asymptotic value of Ta- 
in principle, experiments hint that ^^'^Cs molecules may 
be sufficiently stable against decay caused by inelastic 
collisions 0, 0] but we find that the crossover will never- 
theless be subject to a mechanical instability due to the 
attractive interaction between the molecules. Inelastic 
losses can also be induced in the vicinity of the Feshbach 
resonance by the presence of Efimov states. Indeed, by 
using methods employed in Ref. [1^, we find that these 
three-body states lie in the narrow magnetic field range 
of\B- Bo\<Om G. 

The mechanical instability may be overcome on the 
BEG side by a positive background scattering length for 
the molecules, however it is not possible to properly treat 
this in a single-channel NSR setting and the inclusion 
of this and other effects of molecular interactions is an 
interesting problem for future exploration in the context 
of a two-channel model. 



IV. OUTLOOK AND CONCLUSIONS 

The effects described in this paper may be readily in- 
vestigated in a variety of ways. Pairs of bosonic atoms 
can be experimentally detected using techniques that 
have been successfully employed to observe fermionic- 
atom pairs, for example, by the measurement of cor- 
relations in atom shot noise [2^ or by radio-frequency 
spectroscopy [23j . Furthermore, the experimental obser- 
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vation of a pronounced reduction in the critical temper- 
ature for condensation of paired and unpaired atoms is 
in itself compelling evidence that pairing is taking place. 
As is shown in Figs. \5\ and [51 this striking effect can in- 
deed be seen at relatively high temperatures where the 
gas remains mechanically stable. 

Although the mechanical instability of the cesium gas 
may be overcome on the BEC side of the resonance due to 
a positive background scattering length for the molecules 
as mentioned in the previous section, it remains to be 
shown experimentally whether the gas is stable enough 
against inelastic losses. Indeed, measurements performed 
in Ref. [Q] found inelastic losses due to dimer-dimer col- 
lisions of 1 X lO^^^cni'^/s around i? = 46 G which then 
increase for a fixed temperature as Bq is approached. On 
the other hand, the loss rate was found to be strongly 
temperature-dependent and this may serve to signifi- 
cantly stabilize the gas near Tc. 

Although not the focus of this paper, the condensa- 
tion of unpaired atoms can also be used to experimen- 
tally shed light on some other interesting physics. For 
example, below Tc a transition between bound and un- 
bound (half-)vortex pairs will take place as the line is 
crossed from the PC -I- AC phase to the PC phase [r3.[l3j. 
It would be very exciting if this could be directly observed 
as in the case of the Kosterlitz-Thouless transition in a 



two-dimensional Bose gas [2J1. However, it is not possible 
to make a definite prediction regarding the mechanical 
stability of the gas across this phase transition within the 
framework presented here [25]. Further insight may be 
obtained by using a two-channel approach to treat molec- 
ular interactions beyond mean-field theory. This has the 
added advantage that a background scattering length for 
the molecules as well as atom-molecule interactions can 
be readily included, and would be an interesting topic for 
further research. 
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